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Abstract 

We present a fully self-consistent and relativistic calculation of the nucleon 
mean free path in nuclear matter and finite nuclei. Starting from the Bonn 
potential, the Dirac-Brueckner-Hartree-Fock results for nuclear matter are 
parametrized in terms of an effective a-co Lagrangian suitable for the rela- 
tivistic density-dependent Hartree-Fock (RDHF) approximation. The nucleon 
mean free path in nuclear matter is derived from this effective Lagrangian tak- 
ing diagrams up to fourth-order into account. For the nucleon mean free path 
in finite nuclei, we make use of the density determined by the RDHF calcula- 
tion in the local density approximation. Our microscopic results are in good 
agreement with the empirical data and predictions by Dirac phenomenology. 
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I. INTRODUCTION 



One important quantity of a medium is the mean free path of its elementary constituents. 
In nuclear physics, the nucleon mean free path is of special importance since it can be large 
compared to the nuclear size such that the basic assumption of the independent particle 
motion of the shell model is reasonable |I| • In nuclear reactions, the nucleon mean free path 
is a useful concept for summarizing a large number of experimental data [0,0. Furthermore, 
in the investigation of heavy-ion reactions, the mean free path of nucleons and other hadrons 
are often used to estimate the number of two-body collisions and the reabsorption effect of 
the medium on the production cross sections of hadrons |§ . 

Most calculations of the nucleon mean free path have been done in the framework of 
non-relativistic dynamics |@-11|, based on e.g. the phenomenological Skyrme force. Char- 
acteristic for the early theoretical investigations ||[J is the underestimation of the nucleon 
mean free path by up to a factor of two as compared to the empirical value. The proper 
treatment of the nonlocality of the nucleon optical potential resolves much of the discrep- 
ancy between the theoretical prediction and the empirical data |7]||. The nonlocality of 
the nucleon optical potential (or mean field) leads to the reduction of the nucleon mass and 
consequently increases the nucleon mean free path, which is known as the Negele-Yazaki 
enhancement j7|. 

Recently there have been some relativistic calculations of the nucleon mean free path 
fl2HT5l, based on either Dirac phenomenology or the relativistic impulse approximation for 
the nucleon optical potential. In this paper, we present a fully self-consistent and relativistic 
calculation of the nucleon mean free path in nuclear matter, starting from the Bonn potential 
as the realistic nucleon-nucleon (NN) interaction. The nucleon self-energy (optical potential) 
is derived from the Dirac-Brueckner-Hartree-Fock (DBHF) results for nuclear matter which 
includes the important medium effects. 



This work is a continuation of our effort fl6HT9| to describe self-consistently the prop- 
erties of nuclear matter, finite nuclei and nuclear reactions based on the same realistic NN 
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interaction. There are two aspects to this problem. First, one needs a realistic NN interac- 
tion which is ultimately determined by the underlying dynamics of quarks and gluons and 
should in principle be derived from quantum chromodynamics (QCD). However, due to the 
nonperturbative character of QCD in the low-energy regime relevant for nuclear physics, we 
are far away from a quantitative understanding of the NN interaction in this way. On the 
other hand, there is a good chance that conventional hadrons, like nucleons and mesons, 
remain the relevant degrees of freedom for a wide range of nuclear physics phenomena. In 
that case, the overwhelming part of the NN potential can be constructed in terms of meson- 
baryon interactions. In fact, the only quantitative NN interactions available up till now are 
based on meson-exchange; a well-known example is the Bonn potential |16| , |20[| which we 
apply in this work. 

The second aspect of the problem concerns a suitable many-body theory that is able to 
deal with the bare NN interaction which has a strong repulsive core. The Brueckner approach 
2T^p3[ and the variational method p^ , p5| have been developed for this purpose. However, 



when using two-body forces, both many-body theories are not able to reproduce correctly 
the saturation properties of nuclear matter. Inspired by the success of Dirac phenomenology 
in intermediate-energy proton-nucleus scattering |2?||27] and the Walecka model (QHD) for 
dense nuclear matter P8|j29 |, a relativistic extension of the Brueckner approach has been 



initiated by Shakin and co-workers |3(| , frequently called the Dirac-Brueckner-Hartree-Fock 



(DBHF) approach. This approach has been further developed by Brockmann and Machleidt 
|16| , p!7[| and by ter Haar and Malfliet |31[j. The common feature of all DBHF results is that a 



repulsive relativistic many-body effect is obtained which is strongly density dependent such 
that the empirical nuclear matter saturation can be explained. The Bonn potential and the 
DBHF approach thus provide a reasonable starting point for pursuing the longstanding goal 
of self-consistently describing nuclear matter, finite nuclei and nuclear reactions based on 
the same realistic NN interaction. 

In order to carry out a systematic self-consistent study of nuclear properties, one usually 
parametrizes the DBHF results for nuclear matter in terms of an effective Lagrangian which, 



in relativistic density-dependent Hartree Fock (RDHF) approximation, leads to the same 
predictions for nuclear matter as the original DBHF calculation. The effective Lagrangian, 
with its parameters determined by the underlying NN interaction, can then be used in other 
domains of nuclear physics, e.g. the structure of finite nuclei and nuclear reactions. Different 
schemes for this parametrization have been proposed [P^-pB| . We use in the present work the 



scheme recently suggested by Brockmann and Toki [[34]] in which the DBHF results for nuclear 
matter with the Bonn potential are parameterized in terms of an effective a-u Lagrangian. 
This scheme is however extended to the relativistic density-dependent Hartree-Fock (RDHF) 
approximation, which is more appropriate than the relativistic density-dependent Hartree 



(RDH) approximation of Ref. [34]. The coupling constants of these effective mesons are 
density dependent, and are determined from the underlying bare NN interaction via a DBHF 
calculation. The nucleon self-energy (optical potential), and hence the nucleon mean free 
path, are then calculated based on this effective Lagrangian up to the fourth order Feynman 
diagrams. For the calculation of the nucleon self-energy and mean free path in finite nuclei, 
we use the nucleon density determined by the RDHF calculation [3B|] in the local density 
approximation. 

We outline the formalism of this work in section II. The results and discussion are 
presented in section III. The paper ends with a brief summary in section IV. 

II. FORMALISM 

The relativistic Bonn potential to be used in this work is constructed in terms of the 
Thompson equation which is a three-dimensional reduction to the original Bethe-Salpeter 
equation. The kernel of the Thompson equation, V(q', q), is the sum of one- meson-exchange 



amplitudes of certain bosons with given mass and coupling. In the OBE Bonn model |16 
six nonstrange bosons with mass below 1 GeV are used. Three sets of potential parameters, 
denoted by Bonn A, B and C, have been proposed (and are given in Table A. 2 of Ref. 
|nj ) . The main difference between the three parameter sets is the cut-off mass for the ttNN 
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vertex, which is 1.05, 1.2 and 1.3 GeV for Bonn A, B and C, respectively. Consequently, 
the three potentials differ in the strength of their tensor force component; Bonn A has the 
weakest tensor force. All three potentials reproduce the deuteron properties and the phase 
shifts of NN scattering accurately (cf. Refs. fl6| , [lT|j ). 

The Bonn potential is used in the DBHF calculation for nuclear matter. The essential 
point of this approach is the use of the Dirac equation for the description of the single-particle 
motion in the nuclear medium 

[a k + /3(m + U s ) + U v ]u(k, s) = Eu(k, s) (1) 

where Us is an attractive scalar field and U v is (the timelike component of) a repulsive 
vector field. 

As in conventional Brueckner theory, the basic quantity in the DBHF approach is the G- 
matrix which satisfies the in-medium Thompson equation (also known as relativistic Bethe- 
Goldstone equation) |i6||i7| , |T9 



G(q',q|P,*) = VV,q)+/ ^Lv^)*?— 9S^ <2(k, q|P, z) (2) 

J \ Z7V ) £/ (l/2)P+k Z ^(l/2)P+q - ^(l/2)P+k 

with 

fh = m + U s and £ k = (k 2 + m 2 ) (1/2) 

and P the cm. momentum of the two colliding nucleons in the nuclear medium. 

Since the kernel of the in-medium Thompson equation depends on the solution of the 
Dirac equation, while the Dirac equation needs the scalar and vector potentials which are 
determined from the G-matrix, one is dealing with a self-consistency problem fll6yi7f . The 
nuclear matter properties are then obtained from the in-medium two-body interaction, the 
G-matrix. Applying the Bonn A potential, the DBHF calculation predicts that nuclear 
matter saturates at 0.185 fm -3 with an energy per nucleon S/A=-15.Q MeV, which is in 
good agreement with the empirical values. More results and discussions concerning the 
properties of nuclear matter as predicted by the DBHF approach can be found in Refs. 



As proposed by Brockmann and Toki [34], the DBHF results for nuclear matter can be 



parametrized by an effective Lagrangian, in analogy to the a-u model of Walecka: 

£ = #7^" - m - g a (p)(f>a - gMlM^ 

+ \{d^?-\m 2 Al - \{d^ - d v & + \ml^ (3) 

where ip is the nucleon field, while cj) a and (p^ are the effective sigma and omega fields, respec- 
tively. The masses of the effective sigma and omega mesons are kept fixed at their values in 
free-space scattering. However, the density-dependent coupling constants are choosen such 
as to reproduce the DBHF results for nuclear matter when Eq. (3) is applied in the RDHF 
approximation. 

Treating the effective coupling constants locally as numbers and calculating in the RDHF 
approximation, the nucleon self-energy can be expressed as 

E(^) = £ 5 (M + 7°£o(M + 7"kXy(M (4) 

where £5, S and Ey denotes the scalar component, the time-like part of the vector compo- 
nent and the space-like part of the vector component of the nucleon self-energy, respectively. 
Explicitly the real part of the nucleon self-energy is given by [[29| , j3~7| 

2 g 2 M f k - 



2 a (p) r k F r,m„ 
E s (p,k) = --^ dqq^ 
7T mi Jo hi* 



+ Wkio d ^^9l(p)eAk,q)-^(p)e^k, q )] (5) 

W m 2 9 2 M f kF . 2 

So(P,«) = -o 5~ / dqq 

7r z m^j Jo 

+ tAi: I"" dqq[gl(p)QAk, q) + 2&(p)e u (k, q)] (6) 

lbTf^K Jo 

E v(P> k ) = -g^a dqq^[9l(p)^ ?) + 2fl£(p)*«(*, q)\ (7) 



where 



<S>i(k,q) = In 



Mk,q) + 2kq 
Ai(k,q)-2kq 



®i(k,q) 



A t (k,q)Q i (k,q) 
Akq 



- 1 



Ai{k, q) = k 2 + q 2 + m 2 - (q - k ) 2 , 



i = 



and 



m* k = m + E s (p, k) 



(8) 



k* = k(l + E v (p, k)) 



(9) 



k = (k* 2 + m* 2 ) 1 / 2 + S, 



(10) 



In determining the effective coupling constants from these expressions, we drop the space- 
like component of the vector potential and calculate at the Fermi surface, by identifying Eg 
and So with Us and If® obtained in the DBHF calcaultion. This is a reasonable assump- 
tion since the space-like component is very small and the potentials are only very weakly 
momentum dependent. 

There are mainly two differences between our effective Lagrangian, Eq. (3), and the 
Walecka model |28U29"f . First, the coupling constants in our model are no longer free param- 
eters fitted to the nuclear matter saturation properties; these effective coupling constants are 
determined by the DBHF calculation in which a realistic NN interaction is used. Second, the 
coupling constants of our effective Lagrangian are density dependent, whereas those in the 
Walecka model are density independent. The absence of density dependence in the Walecka 
model may be responsible for its unrealistically large incompressibility. 

In the RDHF approximation, the real part of the nucleon self-energy contains the energy- 
independent Hartree contributions (see Fig. la) as well as the energy-dependent Fock con- 
tributions (see Fig. lb). The lowest order contribution to the imaginary part of the nucleon 
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self-energy is the fourth-order Feynman diagram which is characterized by two-particle-one- 
hole (2plh) intermediate states (see Fig. lc). The nucleon lines in these Feynman diagrams 
are described by dressed nucleon propagators, which corresponds to performing the calcu- 
lation on the Hartree-Fock ground state and taking account of all Hartree-Fock insertions. 
The explicit expressions for the imaginary part of the nucleon self-energy have been given 
in Ref. ||38|| . The derivation of the nucleon self-energy from the Walecka model has been 



discussed in detail in Refs. [p^ , p8[| . For the effective Lagrangian used in the present work, 
the expressions for the nucleon self-energy are the same, but with the coupling constants for 
sigma- and omega-exchange replaced by the density-dependent ones as determined in the 
nuclear matter DBHF calculation. 

The space-like part of the vector potential, Ey, is rather small compared to other terms 
in eq. and can be absorbed into the scalar potential and the time-like part of the vector 
potential by the following transformation 

Us= 1 + Ey ' Uv= 1 + (U) 
In terms of the scalar potential, Us, and the vector potential, Uy, the momentum of a 
nucleon propagating through a uniform nuclear medium can be determined from 

E= [{m + Usf + k^ + Uv (12) 

This can be rewritten as 

^ + V + iW = E-m+ {E ~ m)2 (13) 
2m 2m 

with 

V = U SR + U VR + ^ E ~ m) U V R + ^-(UIr + U VI - U 2 SI - Uy R ) 

m 2m 

W = U SI + U VI + ( E ~ m) U VI + -{UsrUsi - U VR U VI ) 

m m 

where we distinguish between the real and imaginary part of the scalar and vector potential 
given by 

U s = U SR + iU S i 
8 



U v = U V r + iUyi 

Since (V+iW) can be identified as the Schrodinger equivalent potential which is the nucleon 
optical potential in the non-relativistic approach, Eq. fll3|) is identical to the non-relativistic 
dispersion relation, except for the relativistic correction, (E — m) 2 /2m. 

Since the potentials Us and Uy are complex, the nucleon momentum is also complex and 
can be expressed as 

k = k R + ikj (14) 

The nucleon mean free path, A, is related to the imaginary part of the nucleon momentum 
by0 

A - W, < 15 > 



From Eqs. (|i3|), dl4|) and (jig) we obtain an analytical expression for the nucleon mean 
free path 

A = -{-m(E — m — V+ {E - ^ - ) + m[(E — m-V+ ^ E ~ m ^ f + w*ii/ 2} -i/2 (16) 
2 2m 2m 

Since empirical information on the nucleon mean free path is usually obtained by 
analysing nucleon-nucleus scattering data, it is also of interest to perform a microscopic 
calculation of the nucleon mean free path in finite nuclei. To calculate the nucleon self- 
energy (optical potential) and mean free path in finite nuclei, one often makes use of the 
local density approximation. With this approximation, the spatial dependence of the nu- 
cleon self-energy is directly related to the density of the nucleus under consideration. Thus 
in addition to the expressions for the nucleon self-energy in nuclear matter, as outlined 
above, we also need to know the density of the finite nucleus, so that we can calculate the 
nucleon mean free path in finite nuclei. To attain self-consistency of our calculations, the 
density of the finite nucleus must be determined in a RDHF calculation with the effective 



Lagrangian, eq. (3). Such a calculation has recently been carried out by Fritz et al. |36 
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for 40 Ca, which we use in our calculation of the nucleon mean free path. Note that in finite 
nuclei, the incident energy TJ ab is related to the total energy E by 

F _ m 2 + m T (m + T lab ) 

~ ((m + m T y + 2m T T lab y/i [ j 

where m and m T are the masses of the nucleon and the nucleus, respectively. 



III. RESULTS AND DISCUSSIONS 

The coupling constants of the sigma and omega meson of the effective Lagrangian are 
shown in Fig. 2a (sigma) and 2b (omega). Both effective coupling constants drop with 
increasing density. There are some differences between the effective coupling constants 
derived from the Bonn A, B and C potential; those based on Bonn C decrease more with 
increasing density. This difference can be traced back to differences in the tensor-force 
strength of these potentials p9| . 

From the effective Lagrangian we derive the nucleon self-energy (optical potential) up 
to the fourth order Feynman diagrams (see Fig. 1). In Fig. 3, we show the scalar and 
vector potential, Us and Uy, as defined by eqs. (|J) and (|TT|) . We consider two cases with 
density p=(l/2)p (solid curves) and p (dashed curves); p = 0.17 fm~ 3 is the density of 
normal nuclear matter. The results are obtained with the Bonn A potential. The incident 
energy T Zab is related to the total energy E by E = T lab + m. Whereas the real part of these 
potentials depends only weakly on the incident energy and decreases slightly (in magnitude) 
with increasing energy, the imaginary part depends strongly on the incident energy and 
increases (in magnitude) very fast with energy. 

With the scalar potential Us and the vector potential Uy derived from the effective 
Lagrangian, we can calculate self-consistently the nucleon mean free path in nuclear matter 
starting from the bare NN interaction. No parameters beyond those of the Bonn potential 
are involved in the present calculation. The nucleon mean free path is calculated from eq. 
(|IT)|). The results in normal nuclear matter, which simulates the interior of heavy nuclei, 
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are shown in Fig. 4. The solid, long-dashed and short-dashed curves represent the self- 
consistent results based on the Bonn A, B and C potential, respectively, whereas the dotted 
curve is the result derived from the original a-u model of Walecka (QHD-I) ||29|| . The dots 
with error bars represent early empirical values of the nucleon mean free path extracted 
from experiment ||4U|| . The shaded area indicates the estimation of the nucleon mean free 
path based on total reaction cross sections ||. The solid squares correspond to the recent 
empirical data based on Dirac phenomenology in an energy- dependent analysis |l^JT5|] . It 
is clear that there are still large ambiguities in the empirical determination of the nucleon 
mean free path. Within the error range involved, our results based on the Bonn potential 
are in reasonable agreement with the empirical data in the whole energy region considered 
in this work. Comparing with the latest empirical data based on Dirac phenomenology, 
our results are slight larger at low energies and in agreement at high energies, whereas the 
results based on the QHD-I model are in agreement at low energies and too small at higher 
energies. Between the self-consistent results corresponding to Bonn A, B and C, there is 
also some difference which is due to the differences in the effective coupling constants (cf. 
Fig. 2). 

Also of interest is the density dependence of the nucleon mean free path. This information 
is sometimes needed in the description of heavy-ion reactions where dense nuclear matter 
with density up to 2-3po is formed. In Fig. 5 we show the density dependence of the nucleon 
mean free path, corresponding to two incident energies. The results are obtained with the 
Bonn A potential. The mean free path decreases with the increasing density, especially at 
low density. 

Finally we shown in Fig. 6 the nucleon mean free path in 40 Ca as a function of the radial 
distance r. The solid, dashed and dotted curves correspond to nucleon energies T/ ab =150, 
300 and 450 MeV, respectively. In the center of the nucleus, the nucleon mean free path is 
about 2-5 fm, depending on the energy of nucleon. The nucleon mean free path increases 
rapidly at the surface of the nucleus due to the decrease of the nucleon density. 
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IV. SUMMARY AND OUTLOOK 



In this paper, we have presented a fully self-consistent calculation of the nucleon mean free 
path in nuclear matter and finite nuclei, starting from a realistic NN interaction (the Bonn 
potential). In order to facilitate a systematic investigation of nuclear properties, we have 
parametrized the DBHF results for nuclear matter in terms of an effective a-u Lagrangian. 
The density dependent coupling constants of this effective Lagrangian are determined such 
as to reproduce the DBHF results for nuclear matter in the RDHF approximation. The 
effective coupling constants decrease with increasing density. 

In the derivation of the nucleon self-energy (optical potential) based on the effective 
Lagrangian, diagrams up to the fourth order are included. The nucleon mean free path is 
then calculated from the nucleon optical potential through a dispersion relation. Our results 
for the nucleon mean free path in nuclear matter are in good agreement with empirical data. 
We have also analysed the density and energy dependence of the nucleon mean free path 
and find that it decreases with increasing density and energy. 

We have also calculated the nucleon mean free path in finite nuclei by means of the local 
density approximation. For this purpose, we made use of the density of 40 Ca as determined 
in a RDHF calculation with the effective Lagrangian, eq. (3) [jSEfl . In the center of the 
nucleus, the nucleon mean free path was found to be in reasonable agreement with empirical 
data. 

Another aspect is the temperature dependence of the nucleon mean free path, in addition 
to its dependence on density and energy. This is of interest for heavy-ion reactions where a 
piece of hot and dense nuclear matter is formed. This problem is under investigation. 
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FIGURES 

FIG. 1. Feynman diagrams for the calculation of the nucleon self-energy in nuclear matter, (a) 
Hartree diagram, (b) Fock diagram, (c) fourth order diagram. 

FIG. 2. Density dependence of the coupling constants for (a) sigma and (b) omega. 

FIG. 3. Energy and density dependence of the scalar potential Us ((a) and (b)) and vector 
potential Uy ((c) and (d)) in nuclear matter, using the Bonn A potential 

FIG. 4. Energy dependence of the nucleon mean free path in normal nuclear matter. The solid 
dots with error bars are from Ref. [40] , the shaded area is from Ref . [3] and the solid squares from 
Ref. [14] 

FIG. 5. Density dependence of the nucleon mean free path at two energies. Results are obtained 
with the Bonn A potential 

FIG. 6. Nucleon mean free path in 40 Ca as function of radial distance r. The solid, dashed and 
dotted curves correspond to 7) a 6=150, 300 and 450 MeV, respectively. The results are obtained 
with the Bonn A potential. 
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